Instability driven fragmentation of nanoscale fractal islands 
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Formation and evolution of fragmentation instabilities in 
fractal islands, obtained by deposition of silver clusters on 
graphite, are studied. The fragmentation dynamics and sub- 
sequent relaxation to the equilibrium shapes are controlled by 
the deposition conditions and cluster composition. Sharing 
common features with other materials' breakup phenomena, 
the fragmentation instability is governed by the length-to- 
width ratio of the fractal arms. 

PACS numbers: 68.55.Jk, 61.43.Hv, 81.16.Rf, 68.55.-a 



The formation, evolution and consequences of instabil- 
ities in materials systems have been the subject of fruit- 
ful inquiry for a long time |l| [I| . These investigations 
have been performed in various disciplines (e.g. hydro- 
dynamics, solid state physics and materials science) and 
have been applied to systems characterized by different 
length-scales, from the macroscopic |l]-|4j to the nanoscale 
H regimes. Such studies provide deep insights into fun- 
damental issues pertaining to the interplay and balance 
governing the bulk and surface contributions to the ener- 
getics of materials, their morphological stability, and the 
relaxation dynamics in non-equilibrium complex systems. 
Furthermore, studies of such materials instabilities have 
important technological implications concerning the for- 
mation and control of materials systems of desired prop- 
erties (e.g. nanojets ||, nanowires ||, thin films and 
other surface supported structures). 

Here we report on experimental and theoretical inves- 
tigations of the development and evolution of morpho- 
logical instabilities in supported fractal films, generated 
through the deposition of a narrow size-distribution of sil- 
ver clusters on a graphite substrate at room temperature 
We find that such deposited structures fragment 
(at ambient temperature) through a surface self-diffusion 
mechanism [4b], and consequently the mass density and 
its fractal dimension are conserved throughout the pro- 
cess. We demonstrate that the fragmentation dynam- 
ics and the relaxation to the equilibrium shape, can be 
controlled by the deposition conditions as well as by the 
composition of the incident clusters. Analysis of the fully 
fragmented fractal islands leads to identification of the 
ratio between the length and the width of the fractal 
branch (prior to fragmentation) as the critical parameter 
characterizing the instability that underlies the fragmen- 
tation. The emergence of such a characteristic parame- 
ter as well as its critical value (typically of the order of 



4.5), implies that the fragmentation of surface-supported 
solid fractal structures studied here shares the same uni- 
versality class as other instability driven fragmentation 
phenomena 

Under appropriate conditions, films grown through 
the deposition of preformed clusters exhibit a non- 
equilibrium fractal pattern [[7]||. The fractal shape of 
the supported islands is a consequence of the diffusion 
limited aggregation (DLA) mechanism JToj] , and studies 
pertaining to the morphological stability of the films pro- 
vide valuable information about the intrinsic properties 
of these non-equilibrium objects. However, experimental 
reports about fractal stability are scarce |TT|JT^ ] . In these 
experiments, fractal islands relax upon thermal anneal- 
ing into more compact forms and no fragmentation has 
been observed. Nevertheless, numerical simulations pre- 
dicted a late (post-growth) fragmentation stage, with pe- 
riphery diffusion as the dominant mass transport mech- 
anism Our experiments provide evidence for 
fractal fragmentation within a time-scale that depends 
strongly on the fractal branch width and on the surface 
self-diffusion coefficient. The surface diffusion current is 
proportional to the local gradient of the fractal branch 
curvature ( [gjUJl^], and Eq. (||) below). 

Gas-phase neutral silver clusters, produced by a gas- 
aggregation cluster source Q, are deposited at low im- 
pact energy (0.05 eV/atom), on a room temperature 
cleaved graphite surface maintained at a high vacuum 
(10~ 9 torr). The neutral cluster size-distribution is mea- 
sured by a time-of-flight mass spectrometer. The distri- 
bution is characterized by a Gaussian with a peak at 150 
atoms (i.e. clusters with a mean diameter of 2 nm) and a 
full width of half-maximum corresponding typically to 20 
% of the mean diameter value. A crystal quartz micro- 
balance measures the flux of the incident neutral clusters 
(of the order of 10 10 clusters/cm 2 s, with deposition times 
between 1 to 30 minutes) . Since the per-atom kinetic en- 
ergy of the incident clusters is very low compared with 
their binding energy (1.2 eV/atom) the clusters migrate 
on the surface as a whole and grow into islands j?J which 
can continue to evolve past the formation stage. The pre- 
pared samples were either transferred in air and imaged 
in a scanning electron microscope (SEM), or transferred 
into high vacuum and imaged with a non-contact Atomic 
Force Microscope (AFM). The similarity between the im- 
ages obtained in both cases, shows that the transfer in air 
does not affect the island morphologies. We verified by 
X-ray photoemission spectroscopy that further evolution 
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FIG. 1. SEM images of silver fractal islands and his- 
tograms of their branch- lengths (on the right), obtained from 
deposition of Agiso clusters on graphite, for two different cov- 
erages: (a) 6 ML, and (b) 2 ML. 

of the deposited islands after transfer in air is hindered by 
surface oxidation. Therefore, we can define an observa- 
tion time window (OTW) that starts from the beginning 
of the deposition and ends when the sample is exposed to 
air (typical observation time is of the order of one hour). 

In previous studies we have shown that the morpholo- 
gies of islands generated by cluster deposition are the con- 
sequence of the interplay between the coalescence time 
(that is, the time for an added cluster to merge into the 
island) and the aggregation time interval, t a , which is the 
time-lapse between two successive arrivals of clusters to 
the island M. When the coalescence time exceeds t a a 
ramified fractal shape develops with a constant branch 
width Iq. This value depends on the size of the inci- 
dent clusters, on the flux and, to some extent, on the 
deposition time at high coverage, lo decreases for larger 
clusters and/or higher flux, and it increases with time at 
high coverage 0. 

In the first set of experiments we illustrate the effect 
of the fractal branch-width on the fragmentation propen- 
sity. Figs.|l]aand[|b show SEM images of typical individ- 
ual Ag fractal islands (with the fractal dimension |l(]] 
Df = 1.7). The two coverages that we show were ob- 
tained by carrying out the deposition for 7 minutes, and 
subsequently blocking half of the substrate by a mask and 
continuing the deposition onto the unmasked half for an 
additional interval of 20 minutes. The islands shown in 
Figs [l|a and [l]b have approximately the same global size 
(L ~ 1 fim). At large coverage (6ML, unmasked part), 
the grown islands are characterized by thick branches 
(lo = 25 nm, see Fig. |l]a) and they are not fragmented, 
except occasionally at the center of the island. On the 
other hand, at small coverage (2ML, masked part) the 
fractal islands show thinner branches (Iq = 15 nm, see 




FIG. 2. Influence of an impurity in the incident clusters 
on the stability of the fractal islands. The SEM image in (a) 
corresponds to deposition of Agioo clusters and those in (b) 
and (c) correspond to deposition of such clusters but with 
0.5% and 1% oxygen impurity, respectively. In all the images 
a Voronoi polygon construction is superimposed. On the right 
of the images analysis[19] of the fractal dimension is displayed, 
and the asymptotic values are given. 

Fig. [j](b)) and, within our OTW, they exhibit fragmen- 
tation throughout the island. Treating the fractal as a 
set of consecutive bars and focusing on the main branches 
(defined as the longest continuous quasi-linear segments), 
the histograms on the side of the images in Fig. [I] de- 
pict the number distributions of bars as a function of 
the bar length, showing a shift to smaller length as the 
branch width decreases. For the system shown in Fig. ^3 
the fragments have not reached their equilibrium spher- 
ical shape. The non-observability of fragmentation in 
Fig. Qa and the incomplete fragmentation in Fig. 0b, 
are consequences of our OTW value and the slow room- 
temperature fragmentation kinetics (see below). 

In a second set of experiments we vary the composition 
of the incident clusters by adding trace amounts of oxy- 
gen and water during the silver cluster formation stage 
(i.e. before deposition). Fig. || shows SEM images of 
(multi-island) fractal films obtained, for the same cover- 
age, through deposition of pure silver clusters (Fig. ^|a), 
and via deposition of clusters with varying trace amount 
of oxygen (Figs. ||b and |^c). An overall characteriza- 
tion of the deposited film can be obtained through a 
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2D Voronoi polygon construction |L8| , with each Voronoi 
cell being associated with a single island. The den- 
sity of Voronoi polygons is determined by the nucleation 
sites (most likely defects) on the graphite surface (10 s 
sites/cm 2 ). As aforementioned, islands obtained from 
deposition of pure clusters under the above conditions 
do not undergo fragmentation (Fig. |^a), whereas those 
resulting from clusters with variable trace amount of oxy- 
gen are fragmented. The measured density of fragments 
is the same for the data shown in Fig. and |]c, and 
in the latter the shape relaxation of the fragments has 
been completed (resulting in spherical shape). Remark- 
ably, rather similar asymptotic values of the mass fractal 
dimension (D f = 1.70 ± 0.10, see right panels 

in Figs. ||(a-c)) were determined for both the pure and 
impurity-containing cluster depositions. 

From the above observation of a common fractal di- 
mension, the growth of a continuous fractal island should 
precede its subsequent fragmentation into many pieces. 
The fragmentation process was explored through exami- 
nation of the correlation function between the fragments 
belonging to a given island. Such an analysis (per- 
formed for numerous islands) reveals a nearly constant 
nearest-neighbor distance A between fragments belong- 
ing to the same fractal branch. The histogram of the 
nearest- neighbor distance peaks at A = 75 ±10 nm. From 
the size distribution of the fragments we deduced a mean 
fragment diameter 2R = 30 ± 3 nm. Treating the frac- 
tal as a set of uniform cylindrical branches, the width Iq 
of a branch that would have existed prior to fragmen- 
tation can be estimated from the conservation of mass, 
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3A/n- The above considerations lead to a 



ratio X/l() = 4.8. This value is in good agreement with 
that obtained from linear stability analysis, yielding an 
optimal wavelength of 4.443Zo [4c] for fragmentation of a 
solid cylinder induced by a capillary instability (through 
a surface self-diffusion mechanism), and interestingly it 
is close to the value predicted by Rayleigh for breakup of 
a cylindrical column of an inviscid liquid [Q. 

To further investigate the growth and the breakup 
mechanism we performed simulations [^of of the (axial) 
growth of a material with a 3D cylindrical symmetry, 
by adding spherical clusters to each side of a cylindrical 
structure (lying along the z axis) at an aggregation time 
interval t a . The evolution of the shape of the axisym- 
metric material structure [4b] is described by 
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. Here dh/dt and s are, respectively, the local transver- 
sal displacement per unit time and the curvilinear coor- 
dinate of the contour of the material structure on the 
p — z plane of the cylindrical coordinate system (p is the 
local radius of the cylinder along the z axis), k is the 
local curvature expressed as k = K\ + K2, where K\ and 
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FIG. 3. Results from 3D simulations, whose procedure is 
described in the text. In the bottom panel the shape of the 
evolving cylindrical structure and fragmentation events are 
depicted versus time. In the upper panel the length of the 
longest continuous part of the cylindrical structure (i.e. the 
intact length) is plotted versus the simulation time; each frag- 
mentation event is signaled by the appearance of a saw-tooth 
shaped discontinuity, and arrows indicate a change in the clus- 
ters deposition flux. The addition of clusters is stopped for 
t > 78 tu (rightmost arrow) The length and time units are 
scaled as described in the text. 



«2 are the curvatures defined in the p — z plane and in 
the plane normal to that which is tangent to the surface 
of the material structure, respectively. The fourth-order 
diffusion constant B [Q is given by B = D s jfl 2 So/kBT 
where D s is the surface self-diffusion coefficient, 7 is the 
surface tension, f2 is the atomic volume, and So the num- 
ber of atoms per unit surface area. The length and time 
units (tu) are selected by setting both the diameter of the 
incident spherical clusters, and the value of B, to unity. 

For a constant aggregation time interval t a an elon- 
gated cylinder of nearly constant width, with semispher- 
ical caps at both ends, develops at the early growth stage. 
When the length of the structure exceeds a critical value 
an instability develops, resulting in fragmentation into a 
string of equidistant spheres (in agreement with the 2D 
simulations in ref. fl4|| ). 

To model low-frequency experimental fluctuations of 
the cluster source, the flux of the incident clusters is var- 
ied in the simulations during growth. Typical results of 
such a simulation are shown in Fig. ^ In this particular 
simulation we used an aggregation time t a =0.1 tu for 
an initial period of 22 tu, followed by a period of 5 tu 
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during which t a was reduced to a value of 0.09 tu, and 
subsequently t a was restored to its original value. At the 
early stage (t < 20 tu), an elongated capped cylindrical 
shape of nearly constant width develops. As a result of 
the above flux fluctuation the cylinder develops a local 
perturbation of its width which culminates in fragmenta- 
tion at t — 72 tu. At t — 76 tu we stopped the addition 
of clusters. Subsequent evolution involves the develop- 
ment of morphological instabilities that lead to complete 
fragmentation into a string of spheres, which are not rig- 
orously equidistant. It is of interest to remark on the evo- 
lution of the three fragments for t > 90 tu. We observe 
that while eventually all three underwent binary frag- 
mentation, that of the middle one is mass-symmentric 
and those of the two end ones are asymmetric. In fact, 
the morphological instability of the two end-most pieces 
(see the peanut shapes at t = 110 tu) did not achieve 
the critical length-to-width ratio and the process culmi- 
nated in relaxation into a larger-in-size spherical shape. 
Repeating the simulation for several arbitrary tempo- 
ral fluctuations, one obtains a distribution of distances 
A between neighboring spherical fragments, peaking at 
a value of 4.5 Zo , where Iq is the mean diameter of the 
cylindrical structure (prior to fragmentation), in agree- 
ment with the experimental histograms. 

Underlying the evolution of a morphological instability 
of an elongated solid 3D material structure is the depen- 
dence of the surface diffusion rate on the local curvature. 
A linear stability analysis [4c] of Eq. (|l|), limited to a 
small axisymmetric disturbance (with wave- vector fc) of 
the cylinder's radius, gives for the fastest growing dis- 
turbance k 2 pQ — 0.5, where po is the radius of the un- 
perturbed uniform cylinder; i.e., an optimal instability 
wavelength A c = 2 1 / 2 27r j o — 4.443/ - This value of X c /l 
is very close to the aforementioned simulated value of 
4.5Zo- We remark that within linear stability analysis, 
such instability does not develop for a 2D system. 

From the measured mean radius of the relaxed frac- 
tal fragments (Fig. ||c), together with mass conserva- 
tion, wc conclude that both the fragmented and unfrag- 
mented fractal islands are 3D in nature, with a wet- 
ting angle in the vicinity of tt/2 [^l|. The disturbance 
which may culminate in fragmentation can be either of 
thermal origin or caused by substrate induced stresses, 
and the relevant time-scale for fragmentation involving 
surface self-diffusion is r = Iq/B. This implies that 
the kinetics of fragmentation is influenced by Iq and/or 
B. Fig. [l] demonstrates the influence of Iq. Using pa- 
rameters for pure silver p2] , p3| we estimate at 300K, 
t(1 = 25nm) = 2300 s and t(Iq = 15nm) = 300 s. 
Consequently for a typical 2000 s OTW used by us, ob- 
servation of room-temperature fragmentation is (kineti- 
cally) unlikely for the thicker branches (Fig. [l]a), while 
it should be readily observed (at least partial fragmen- 
tation) for the thinner ones (Fig. [j]b). Figs. ||d and |c 
depict complete fragmentation within the same OTW, 



showing evidence for a faster fragmentation kinetics (i.e. 
decrease of r) with increasing trace amounts of oxygen 
in the incident clusters. 

In summary, we found via experiments and theoretical 
analysis, that the room-temperature fragmentation insta- 
bility of fractal islands formed through the deposition of 
silver clusters on graphite, is governed by the ratio of the 
length to the width of the fractal arms; such length-to- 
width control parameter is a common feature in many in- 
stability phenomena In particular, the wavelength 
of the fastest growing instability is determined by this 
critical ratio (having a value of ~ 4.5), and the actual 
observability of fragmentation depends on the fragmen- 
tation kinetic, involving surface self-diffusion, and thus 
on the experimental OTW. The fragmentation dynam- 
ics and the relaxation to the equilibrium shapes can be 
controlled by the deposition conditions (including size se- 
lection of the incident clusters) as well as by the compo- 
sition of the incident clusters. Technological implications 
may include issues pertaining to the growth of surface- 
supported elongated structures (nanowires), and the self- 
organization of structures made of supported nanoparti- 
cles with well-defined sizes and intcrparticle spacing. 

This research is supported by CNRS and the US DOE 
(Grant No. FG05-86ER-45234) . Computations were per- 
formed at the GIT Center for Computational Materials 
Science. 



[1] J.W.S. Rayleigh, proceedings of London Math. Soc, 10, 
4 (1879). 

[2] S. Chandrasckhar, Hydrodynamics and Hydrodynamic 
Stability (Oxford Univ., New York/London 1961). 

[3] J.Eggers, Rev. Mod. Phys. bf 69, 865 (1997). 

[4] (a) W.W. Mullins, J. Appl. Phys. bf 28, 333 (1957); (b) 
F. A. Nichols and W. W. Mullins, J. Appl. Phys. bf 36, 
1826 (1965); (c) F. A. Nichols and W.W. Mullins, Trans. 
Met. Soc. AIME, 233, 1840 (1965). 

[5] M. Moseler and U. Landman, Science 289, 1165 (2000). 

[6] Nanowires, Eds. P. Serena and N. Garcia (Kluwer, Dor- 
drecht 1997). 

[7] B. Yoon, et al, Surf. Sci. 443, 76 (1999). 

[8] I. M. Goldby, et al, Appl. Phys. Lett. 69, 2819, (1996). 

[9] P. Jensen, Rev. Mod. Phys. 71, 1695 (1999). 
[10] T. A. Witten and L. M. Sander, Phys. Rev. Lett. 47, 
1400 (1981). 

[11] R. Q. Hwang, et al, Phys. Rev. Lett. 67, 3279 (1991). 

[12] H. Roder, et al, Nature (London) 366, 141 (1993). 

[13] T. Irisawa. et al, Europhys. Lett. 30, 139 (1995). 

[14] R. Thouy, et al, Phys. Rev. B 56, 5321 (1997). 

[15] R. Sempere, et al, Phys. Rev. Lett. 71, 3307, (1993). 

[16] M.Conti, et al, Phys. Rev. Lett. 80, 4693 (1998). 

[17] J.Eggers, Phys Rev. Lett. 80, 2634, (1998). 

[18] P. A. Mulheran, J. A. Blackman, Phil. Mag. Lett. 74, 



4 



4694 (1995). 

[19] The asymptotic mass fractal dimension was measured for 
large R g by using the power law relation between the 
mass M and the gyration radius R g . 
[20] The numerical method that we used is a 3D extension of 
the 2D simulation described in ref. jwj. Moreover, in ref. 
only fragmentation after growth was simulated. 
[21] M. S. McCallum, et al, J. Appl. Phys. 79, 7604 (1996). 
[22] C.R. Stoldt, et al, Phys. Rev. Lett. 81, 2950 (1998). 
[23] K. Akamatsu, et al, J. Phys. Chem., 104, 10168 (2000). 



5 



